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Introduction 


It  is  well  known  that  the  phase  noise  of  a  sinusoidal  oscillator  is  reduced  by 
the  time-delayed  mixing  scheme  shown  in  figure  la.  This  phenomenon,  called 
the  correlation  effect  (ref  1)  is  important  in  the  phase  noise  analysis  of  radar  and 
communications  systems.  It  is  usually  explained  using  FM  modulation  theory, 
(refs  1  through  4)  where  it  is  shown  that  the  reduction  in  the  sideband  power 
at  the  output  relative  to  the  sideband  power  of  the  source  is  given  by  the  well 
known  formula 


K(u)  =  2[l  —  cos(u;to)]  =  4sin2(arro/2)  (1) 

where  To  is  the  time  delay  and  ll>  is  the  radian  frequency  offset  from  the  carrier. 
A  precise  definition  of  the  function  K{uj)  is  given  below.  In  this  report,  the 
correlation  effect  is  analyzed  using  the  theory  of  stochastic  processes.  The 
signal  model  consists  of  a  fixed  frequency  sinusoid  that  is  phase  modulated  by 
a  random-walk  phase  function,  i.e.,  a  random-walk  phase  noise.  Other  signal 
models  are  possible  but  shall  not  be  considered  in  this  report.  It  is  shown  that  an 
analysis  based  on  this  random  signal  model  leads  to  a  similar,  though  somewhat 
more  detailed  formula  than  equation  (1).  In  one  respect,  the  analysis  presented 
here  places  equation  (1)  on  a  more  rigorous  mathematical  foundation.  In  fight 
of  this  analysis,  it  is  fair  to  conclude  that  equation  (1)  is  a  good  approximation 
which  is  generally  valid  except  for  very  small  frequency  offsets. 

The  problem  considered  here  is  illustrated  in  figure  la.  This  is  a  common 
signal  processing  operation  which  can  be  thought  of  as  the  front  end  of  a  co¬ 
herent  CW  or  FMCW  radar  receiver.  The  signal  v(t)  is  a  sinewave  containing 
a  random  phase  function,  i.e., 

u(f)  =  cos[w0f  +  <p{t)  +  (po]  (2) 

where  (p(t )  is  a  random  process.1  The  phase  noise  originating  in  the  source  oscil¬ 
lator  v(t)  propagates  through  the  mixer  and  low-pass  filter  (LP)  and  manifests 
itself  in  the  output  x(t).  The  problem  is  to  compare  the  phase  noise  spectrum 
of  the  output  signal  x(t )  to  the  phase  noise  spectrum  of  the  source  v(t).  The 
term  phase  noise  spectrum  is  defined  in  the  section  Spectral  Purity.  In  order  to 
simplify  the  comparison  of  these  power  spectra,  it  is  helpful  to  frequency-shift 
the  power  spectrum  of  the  source  oscillator  v(t)  from  the  carrier  frequency  o>o 
down  to  D.C.  This  is  accomplished  very  easily  by  means  of  the  circuit  shown 
in  figure  lb.  In  figure  lb,  the  local  oscillator  is  an  ideal  reference  oscillator 

1For  an  FMCW  radar,  an  FM  modulation  term  should  be  included  in  the  argument  of 
the  cosine.  However,  this  modulation  term  does  not  effect  the  analysis  and  is  omitted  for 
simplicity. 
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u(t)  =  cos  wo t  which  has  zero  phase  noise.  Thus,  the  power  spectrum  of  the 
output  y(t)  is  almost  identical  to  the  power  spectrum  of  the  source  v(t),  except 
that  the  entire  spectrum  is  downshifted  to  D.C.  Observe  that  the  power  spec¬ 
trum  of  v(t)  can  be  easily  measured  in  this  case  by  connecting  the  output  y(t) 
to  a  baseband  spectrum  analyzer.  This  is  the  basic  principle  behind  various 
methods  of  phase  noise  measurement  (ref  5). 

The  main  objective  of  this  study  is  to  compare  the  power  spectral  density 
of  the  mixer  output  x(t)  to  the  power  spectral  density  of  the  source  oscillator, 
i.e.,  y(t).  It  will  be  shown  that  there  is  a  reduction,  i.e.,  a  decrease  in  the  phase 
noise  spectrum  of  the  signal  x(t)  relative  to  the  phase  noise  spectrum  of  y(t). 
The  reduction  in  phase  noise  power  is  characterized  by  the  ratio 

<3) 

where  Sx(u>)  and  Sy(ui)  are  the  power  spectral  densities  of  x(t)  and  y(t),  respec¬ 
tively.  The  quantity  K(u>)  will  be  called  the  reduction  factor  or  the  correlation 
factor.  (In  ref  4  it  is  called  the  delay  function).  Using  dB  notation,  the  reduction 
factor  can  be  expressed  in  the  form 


=  101og10  dB.  (4) 

The  main  goal  of  this  report  is  to  calculate  the  reduction  factor  for  the  random 
signal  model  given  by  equation  (2). 

Spectral  Purity 

The  spectral  purity  of  a  sinusoidal  oscillator  describes  how  closely  the  fre¬ 
quency  spectrum  of  the  oscillator  represents  that  of  an  ideal  sinewave.  The 
term  phase  noise  refers  to  the  short  term  (<  1  second)  frequency  fluctuations 
resulting  from  random  phase  modulation  of  the  sinusoidal  signal.  Phase  noise 
is  also  called  FM  noise.  In  addition  to  phase  noise,  an  oscillator  is  also  subject 
to  amplitude  fluctuations  or  AM  noise.  However,  for  an  amplitude  stabilized 
oscillator,  the  AM  noise  is  usually  much  less  than  the  FM  noise. 

The  power  spectrum  of  a  good  oscillator  is  characterized  by  a  narrow  spec¬ 
tral  peak  at  /  =  /o  (and  one  at  — /o).  A  simple  measure  of  spectral  purity  is 
given  by  the  linewidth  of  the  power  spectrum  S(f)  as  shown  in  figure  2.  Note 
that  the  power  spectrum  as  used  here  is  defined  for  both  positive  and  nega¬ 
tive  frequencies.  Moreover,  it  is  an  even  function  of  frequency.  The  oscillator 
linewidth  is  characterized  by  either  the  full-width  at  half-maximum  (FWHM) 
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Figure  2.  The  linewidth  of  the  power  spectrum. 
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or  the  half-width  at  half-maximum  (HWHM).  Thus,  an  elementary  measure  of 
spectral  purity  is  given  by  the  ratio 


„  ^3dB  /E\ 

” =  IT’  (5) 

where  /3d b  is  equal  to  the  HWHM.  The  parameter  77  will  be  called  the  linewidth 
parameter. 

A  much  more  precise  measure  of  spectral  purity  is  given  by  the  quantity  £(/), 
called  script-L-F.  It  seems  that  £(/)  was  first  introduced  into  the  engineering 
literature  by  V.  Van  Duzer  in  1964  (ref  6).  Another  early  appearance  of  C(f) 
in  the  literature  was  due  to  D.J.  Glaze  in  1970  (ref  7).  A  formal  definition  of 
C(f)  was  given  by  D.W.  Allen,  J.H.  Shoaf,  and  D.  Halford  in  1974  (ref  8).  As 
originally  defined  in  reference  8: 

m  =  (6) 

■TO 

where  P(/)  is  the  average  signal  power  in  a  1  Hz  bandwidth  and  at  a  frequency 
offset  /  from  the  carrier,  and  Po  is  the  total  average  signal  power.  To  be  more 
precise,  P(f)  is  defined  to  be  the  power  spectral  density  at  a  frequency  offset  / 
from  the  carrier,  where  /  lies  in  the  range  — /o  <  /  <  00.  Since  it  is  sometimes 
a  point  of  confusion,  it  is  important  to  stress  that  Po  is  not  the  power  spectral 
density  at  /o,  but  the  total  average  signal  power.  It  is  important  to  make  this 
distinction.  An  immediate  consequence  of  the  above  definition  is  that  £(/) 
satisfies  the  normalization  condition 


r  c(t)if =1. 

J-fn 


Therefore,  C(f)  will  be  called  the  normalized  power  spectral  density.  As  one 
might  expect,  £(/)  is  directly  related  to  the  power  spectrum  S(f)  of  the  oscil¬ 
lator.  By  the  definition  of  S(f).  the  total  average  signal  power  is  given  by 


p0  =  r  s(f)df. 

j —00 


Using  the  fact  that  S(f )  is  an  even  function  of  /,  that  is,  S(—f )  =  S(f ),  this 
can  be  written  as  an  integral  over  positive  frequencies  in  the  form 


rOO 

P0=  2  S(f)df. 

Jo 
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The  average  power  contained  in  the  frequency  band  from  f\  to  /2  (0  <  }\  <  f%) 
is  given  by 


Pi2=  [2ZS(f)df 
Jf , 


and  the  fraction  of  the  total  power  which  is  contained  in  that  frequency  band 

Pn_ 

Po  So°°S(f)df  '  1 

The  spectral  density  P(/)  can  be  obtained  from  equation  (9)  or  (10)  by  inspec¬ 
tion:  at  a  frequency  offset  /  from  the  carrier  P(f)  =  2S(fo  +  /).  Therefore, 
one  obtains  the  relation 

c(i)  =  (i2) 

jo 

Consider  now  the  units  (physical  dimensions)  of  £(/).  Because  P(f)  is  a  power 
spectral  density  and  Po  is  a  power,  the  units  of  £(/)  are  Hz-1.  Using  dB 
notation,  £(/)  can  be  written  in  the  form 


£(/)  =  10  log10 


The  units  in  this  case  are  dBc/Hz  which  are  commonly  abbreviated  as  dBc. 

As  defined  by  equation  (6),  the  normalized  power  spectrum  £(/)  is  the  stan¬ 
dard  way  to  characterize  oscillator  spectral  purity.  In  the  measurement  com¬ 
munity,  £(/)  is  called  the  single-sideband  (SSB)  phase  noise.  Measurements  of 
oscillator  spectral  purity  seek  to  measure  £(/)  as  a  function  of  /.  Assuming 
that  the  spectral  profile  is  symmetric  about  fo,  measurements  are  generally 
made  only  for  positive  (or  only  for  negative)  values  of  the  frequency  offset  /. 
For  the  system  designer,  acceptable  limits  of  phase  noise  are  usually  specified 
by  quoting  an  upper  bound  for  £(/)  at  a  specified  offset  /  from  the  carrier. 
For  a  more  detailed  specification,  an  upper  bound  on  £(/)  might  be  specified 
at  several  discrete  frequency  offsets,  /i,/2, •  •  • ,  fn’,  or  over  a  continuous  range 
of  frequencies,  fi  <  f  < 


Stochastic  Phase  Modulation  Theory 

The  analysis  given  here  follows  Papoulis  (ref  9,  2nd  ed.,  p.  321-324).  It  is 
assumed  that  the  signal  source  is  amplitude  stabilized  so  that  the  amplitude 
fluctuations  of  the  oscillator  output  are  negligible.  The  phase  fluctuations  of 
the  oscillator  give  rise  to  a  stochastic  process  of  the  form 

v(t)  =  cos[u>ot  +  <p(t)  +  ipo],  (14) 
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where  the  phase  function  <p(t)  is  a  stochastic  process  and  tpo  is  a  random  variable 
that  is  uniformly  distributed  on  the  interval  [0, 27r].  In  addition,  it  is  assumed 
that  ’-pit)  and  ipo  are  statistically  independent  for  all  t.  There  axe  two  ways 
to  specify  the  process  p(t):  Either  p{t)  is  a  given  stochastic  process  (phase 
modulation)  with  known  statistics,  or  it  is  defined  indirectly  by  the  integral 


¥>(*)=  [  w(s)ds, 
Jo 


where  w(t)  is  a  given  stochastic  process  (frequency  modulation)  with  known 
statistics.  Technically,  equation  (15)  is  a  stochastic  integral  which  must  be 
treated  with  special  mathematical  care;  however,  the  representation  (15)  shall 
not  be  used  here.  For  the  purposes  of  analysis,  it  is  desirable  that  the  oscilla¬ 
tor  signal  v(t)  be  at  least  wide-sense  stationary  so  that  the  process  mean  and 
autocorrelation  function  are  invariant  with  respect  to  a  shift  of  the  time  origin. 

To  simplify  the  mathematics,  the  real  valued  process  (14)  can  be  represented 
as  the  real  part  of  the  complex  process  V (t)  defined  by 

V(t)  =  (16) 

The  statistics  of  v(t)  are  closely  related  to  the  statistics  of  the  process 

$(<)  =  e^+V0\  (17) 

since  V (t)  =  $(f)eIu;ot.  The  objective  now  is  to  find  the  power  spectrum  of  the 
oscillator  signal  v(i).  First  of  all,  the  process  v(t)  has  zero  mean  since 

E\et(f>°~\  =  E[cos<^o]  +  f£^[sint/?o]  =  0, 

and  therefore,  using  the  fact  that  the  <p(t)  and  (p o  are  statistically  independent, 

j  E[V(t)]  =  E(ei^ot+v^)E[eitfio]  =  0 
\£[V*(t)]  =  E{e-i^t+v^)E[e-iv°]=  0, 


which  imply 


£[i>(<)]  =  E 


V(t)  +  V*(f) 
2 


Next,  the  autocorrelation  function  of  the  process  v(t)  can  be  written  in  the  form 

R,(h,t  i)  =  =  e|  | 

=  j{b[ V(t,)V(h)]  +E[V(t2)V*(*1)]  +B[V(<j)V(<0] 

+  E[V*(i>)V(«,)]}.  (18) 


7 


(19) 


Evaluating  the  first  term  yields 

E[V(t2)V(h)]  =  E[eia,o(t2+‘l)ei[^2)+^(<l)]]  £  [e*2*10]  =  0. 

The  second  term  in  equation  (18)  can  be  written  in  the  form 

E[V(t2)V*(h)]  =J^(t2,f1)eiwo(f2-tl),  (20) 

or  equivalently, 

i2v(<2,ii)  =  -R$(t2,<1)eiwo(t2-tl),  (21) 

where 

JM*2,*i)  =  £[#(f2)$*(fi)]  =  (22) 

The  second  and  third  terms  in  equation  (18)  are  obtained  from  the  first  and 
second  terms  by  interchanging  1 1  and  <2 ,  or  equivalently,  by  complex  conjuga¬ 
tion.  Therefore,  substituting  equations  (19)  and  (20)  into  equation  (18),  one 
obtains 

Rv(t2,h)  =  ^[i4(*2,ti)eia'or  +  i24(f2,<i)e-ita'or],  (23) 

where  r  —  t2  — 1\.  This  equation  expresses  the  autocorrelation  function  of  v(t) 
in  terms  of  the  autocorrelation  function  of  $(f).  As  stated  previously,  it  is 
desirable  that  the  random  process  v(t)  be  wide-sense  stationary.  If  R§  is  only  a 
function  of  r  =  then  so  is  Rv.  Thus,  a  sufficient  condition  for  stationarity 

is  easy  to  obtain.  If  the  first  order  p.d.f.  of  the  increment  <^2)— ^(h)  appearing 
in  equation  (22)  depends  only  on  ^  ,  then  R<p  (f 2  •  t\ )  is  only  a  function  of 

r.  Consequently,  equation  (23)  tabes  the  form 

Rv(r)  =  i  [f?$(r)e‘W0T  +  i?S(r)e-*w°r] .  (24) 

Therefore,  the  autocorrelation  function  of  the  process  v(t)  only  depends  on 
the  time  difference  r,  and  since  E[v(t)\  =  0,  this  implies  that  v(t)  is  wide- 
sense  stationary.  In  general,  for  a  Gaussian  phase  modulation,  a  necessary  and 
sufficient  condition  for  the  process  v(t)  to  be  wide-sense  stationary  is  given  in 
appendix  A. 

Having  obtained  the  autocorrelation  function  of  u(t),  it  is  now  possible  to 
obtain  the  power  spectrum  by  simply  taking  the  Fourier  transform  of  Rv(r). 
Using  the  well  known  frequency-shift  theorem,  the  Fourier  transform  of  equation 
(24)  is  given  by 

Sv(w)  =  ^[S$(w-a;o)  +  S$(a;  +  a;o)],  (25) 

where  S<t>(u>)  is  the  power  spectrum  of  $(<).  Of  course,  S^(u>)  is  real  since  by 
equation  (22),  i?$(— r)  =  f?£(r).  Equation  (25)  gives  the  power  spectrum  of 
the  oscillator  signal  u(f)  in  terms  of  the  power  spectrum  of  4>(t).  In  order  to 
proceed  with  the  analysis,  it  is  necessary  to  assume  some  particular  form  for 
the  process  < p{t). 
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Random- Walk  Phase  Modulation 


It  is  known  that  due  to  thermal  and  shot  noise,  the  phase  ip(t)  in  an  electrical 
oscillator  executes  a  continuous  random  walk  (refs  10  and  11).  Therefore,  it 
will  be  assumed  that  (p(t)  is  a  zero  mean  Wiener  process  with  parameter  a. 
The  first  order  p.d.f.  of  tp(t)  is  given  by 


/(¥>>*)  = 


\/2tt at 


D-tp*/2at 


where  t  >  0.  According  to  Edson  (ref  10),  the  coherence  time  rc  is  defined  to 
be  the  time  interval  required  for  the  standard  deviation  (variance)  of  the  phase 
to  increase  by  1  radian.  For  the  Wiener  process  (26),  this  condition  takes  the 
form 

Vax  jV(<)]  =  E[y>2(f)]  =  at  =  1,  (27) 

and  therefore,  the  coherence  time  is  given  by 

tc  =  1/a.  (28) 

One  way  to  define  the  Wiener  process  is  by  the  integral 


(p(t)  =  I  w(s)ds , 

Jo 


where  w(t)  is  a  zero  mean,  stationary,  white  noise  Gaussian  process  with  au¬ 
tocorrelation  function  Rw(r )  =  a5(r).  Thus,  for  a  random-walk  phase  modu¬ 
lation  the  instantaneous  frequency  deviation  u{t)  —  uq  —  w(t)  is  a  zero  mean, 
stationary,  white  noise  Gaussian  process  where  the  rms  frequency  fluctuations 
are  given  by 

<r2  =  Var  [u?(f)  —  a>o]  =  Var  [«>(<)]  =  a.  (29) 

The  variance  is  related  to  the  linewidth  of  the  power  spectrum.  Combining 
equations  (28)  and  (29)  shows  that  the  rms  frequency  fluctuations  are  related 
to  the  coherence  time  by  the  relation 


To  avoid  the  use  of  the  stochastic  integral,  the  Wiener  process  can  also  be 
defined  by  the  following  properties: 


1.  ^(0)  =  0. 

2-  f(<p(t2)  =  <P2\v(ti)  =  <Pi) 


^Kajt-i  - 12) 


(V?2  ~yi)2 
2 a(t2  -h) 


^2  >  ^1  ^  0. 
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3.  The  increments  <p(t2)—<p(ti)  and  <^>(£4)  — <£>(£3)  axe  statistically  independent 
for  every  set  of  points  0  <  ti  <  t2  <  <3  <  t.\. 

4.  The  sample  functions  <p(t)  are  continuous. 

These  properties  have  the  following  simple  consequences: 

5.  The  first  order  p.d.f.  of  the  increment  ip(t  +  r)  —  <p(t)  is  independent  of  t. 

6.  The  joint  p.d.f.  of  <p(ti),cp(t2 ), . . .  ,<p(tn),  is  Gaussian  for  every  sequence 
0  <  h  <  t2  <  •  •  •  <  tn. 

It  is  well  known  (ref  9)  that  the  autocorrelation  function  of  <p(t)  is  given  by 

R<p(t2,h)  —  amin(ti,t2)-  (31) 

This  shows  that  cp(t)  is  not  a  stationary  process.  Nevertheless,  the  process  v(t) 
turns  out  to  be  stationary. 

The  properties  of  the  Wiener  process  will  now  be  used  to  compute  the  power 
spectrum  of  the  oscillator  signal  v(t).  Since  a  linear  combination  of  Gaussian 
random  variables  is  itself  a  Gaussian  random  variable,  the  increment  <^2)  — 
<£>(fi)  has  a  Gaussian  p.d.f.  with  mean  value  zero  and  variance  a|i2  —  <1 1.  Thus, 
the  autocorrelation  function 

J2*(*2,f  1)  =  £(e’M<2)-v(ti)])  (32) 

is  essentially  the  characteristic  function  of  the  Gaussian  random  variable  95(^2)  — 
<p(ti).  From  probability  theory,  if  x  is  a  Gaussian  random  variable  with  mean 
p  and  variance  <r2,  then  the  characteristic  function  of  x  is 

E[eixz]  =  ev*-*2*a/2.  (33) 

Using  this  result  together  with  property  2  of  the  Wiener  process,  it  follows  that 

JU{t2,t1)  =  e-^~t^\  (34) 


or  more  simply, 

^(r)  =  e“  M/2.  (35) 

Taking  the  Fourier  transform  of  this  function  yields  the  power  spectrum 


iS$  (a;) 


(a/2)2  +a>2' 


This  is  the  well  known  Lorentzian  spectrum.  The  half-width  at  half-maximum 
of  the  spectral  peak  is  given  by 


a  1 

"mb  =  2  =  2^: 
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This  relates  the  coherence  time  to  the  the  linewidth  of  the  power  spectrum. 
Using  equations  (36)  and  (25),  the  power  spectrum  of  the  oscillator  signal  v(t) 
is  given  by 


/2)2  +  (u  -  a,,,)2  T  (a/2)2  +  (u  +  u0  f 


This  is  the  desired  result.  The  PSD  (38)  consists  of  two  Lorentzians  symmet¬ 
rically  located  at  —  ±u>o. 


Phase  Noise  of  the  Source  Oscillator 

The  phase  noise  of  the  source  oscillator  can  be  described  by  either  the  RF 
power  spectrum  Sv(u>)  or  by  the  normalized  phase-noise  spectrum  £(/).  These 
two  spectral  representations  are  related  through  equation  (12).  The  average 
power  in  the  signal  v(t),  the  total  signal  power,  is  given  by 

P0  =  E[v2(i)}  =  R„(0)  =  ±  (39) 

where  the  value  1/2  is  obtained  from  equations  (24)  and  (35).  Substituting 
equations  (39)  and  (38)  into  equation  (12),  the  SSB  phase  noise  of  the  source 
oscillator  is  given  by 


£(/)  [(a/2)2  +  u,2  +  (a/2)2  +  (a,  +  2a>„)2  j  ' 


For  a  frequency  offset  u:  such  that  |oj|  <wo,  the  magnitude  of  the  first  term 
is  much  greater  than  the  magnitude  of  the  second  term.  Hence,  it  is  a  good 
approximation  to  simply  write 


4/)  = 


(a/2)2  +l>2J 


Thus,  the  normalized  SSB  phase-noise  spectrum  has  the  same  Lorentzian  spec¬ 
tral  profile  as  the  RF  power  spectrum. 

This  result  can  also  be  obtained  by  means  of  the  mixer  circuit  in  figure  lb. 
In  this  circuit,  the  noisy  source  oscillator  u(t)  is  mixed  with  a  phase-noise-free 
reference  oscillator  u(t)  =  cos(wot)-  The  input  to  the  low-pass  filter  is 


v(t)u(t)  —  cos[wot  +  <p{t)  +  V’o]  cos(u;ot). 


(42) 
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Using  the  cosine  identity 


1  r 

cos  A  cos  B  =  -  [cos(A  —  B)  +  cos( A  +  B)\ ,  (43) 

together  with  the  fact  that  (p(t)  is  a  lowpass  process,  the  filtered  mixer  output 
is  given  by 

y(i)  =  \  cos  [p(t)  +  wo] .  (44) 

Since  ip(t)  is  not  strictly  a  low-pass  process,  this  equality  is  not  exact.  A  random 
process  x(t)  is  a  strictly  low-pass  process  if  =  0  for  |u;|  >  B,  where  B  >  0 

is  the  cutoff  frequency.  Nevertheless,  for  an  oscillator  with  good  spectral  purity, 
the  process  and  therefore  the  process  exp[i<£>(t)],  is  very  nearly  a  low-pass 
process  so  that  equation  (44)  is  a  reasonable  approximation.  To  continue  with 
the  analysis,  let  Y (t)  be  the  complex  signal 

“2 

Note  that  2Y(t)  =  $(t)  as  defined  by  equation  (17).  Proceeding  as  in  section 
2,  it  is  clear  that  B[F(t)]  =  B[F*(t)]  =  0,  and  therefore  E[y{t)\  —  0.  The 
autocorrelation  function  of  y(t)  can  be  written  in  the  form 

=  \{E[Y(ti)Y(h)]  +E[Y(h)Y'(tl)}  +E[Y'(t2)Y(t1)] 

+  B[r*((2)y*((0]}-  (46) 

Using  the  results 

E[Y(h)Y(h)}  =B[F+(<2)r*(t1)]  =0, 

and 

E[Y(t2)Y 'fa)]  =  =  Ij4(t), 

equation  (46)  implies 

Ry(T)  =  (47) 

where  R$  is  defined  by  equation  (22)  and  r  =  t2  —  1 1 .  By  taking  the  Fourier 

transform,  this  immediately  yields  the  power  spectrum 

Sy(“)  =  £$*(«)•  (48) 
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Substituting  the  result  (36)  for  a  random- walk  phase  modulation,  one  obtains 


s»(u)  = 


(a/2)2  +u2\  ' 


This  is  the  desired  power  spectrum.  Equation  (12)  will  now  be  used  to  calculate 
the  normalized  power  spectrum  £(/).  From  equation  (47),  the  total  average 
power  in  the  signal  y(t)  is  found  to  be  Pq  =  1  /8.  Combining  this  with  equation 
(49),  the  normalized  PSD  of  y(t)  is  given  by 


A/)  = 


(a/2)2  +  u;: 


w=27 r/ 


Except  for  an  additional  factor  of  two,  this  is  identical  to  equation  (41).  Conse¬ 
quently,  the  normalized  phase  noise  spectrum  of  y(t)  is  two  times  the  normalized 
phase  noise  spectrum  of  the  source  oscillator  v(t). 

For  the  purpose  of  making  calculations,  the  phase-noise  spectrum  (41)  can 
be  written  in  the  more  convenient  form 


4/)  = 


where  the  HWHM  is  given  by 


7r/3dB  [^  +  (f/f3dn)2 


/3d  B  = 


Alternatively,  this  can  be  expressed  in  dB  notation  in  the  form 


10 logj 


?r/3dB  [1  +  (///3dB)2 


dBc/Hz. 


For  a  36  GHz  oscillator  with  a  linewidth  parameter  rj  =  10~7,  the  normalized 
phase-noise  spectrum  of  v{t),  equation  (53),  is  plotted  in  figure  3.  The  phase 
noise  of  the  signal  y(t)  can  be  obtained  by  simply  adding  3  dB.  The  spectrum 
in  figure  3  is  the  well  known  Lorentzian.  It  is  flat  from  /  =  0  to  /3dB  =  3600 
Hz  and  then  rolls  off  with  a  slope  of  20  dB  per  decade.  In  the  low  frequency 
limit,  as  /  -4  0,  £(/)  -4  -101og(7r/3dB)  =  -40.5  dBc/Hz. 


Autocorrelation  Function  of  the  Mixer  Output 

With  reference  to  figure  la,  the  two  signals  entering  the  mixer  are  the  oscilla¬ 
tor  signal  v(t)  and  the  received  signal  ev(t  —  r0),  where  t0  >  0  is  the  time  delay 
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Figure  3 .  Normalized  power  spectrum  for  a  36  GHz 
oscillator  with  a  HWHM  (3  dB  frequency)  of  3600  Hz 


14 


of  the  received  signal  and  e  >  0  is  its  magnitude.  The  input  to  the  low-pass 
filter  is  the  product 

ev(t)v(t  -  To)  =  £  cosfwof  +  <p{t)  +  ipo]  COs[u>0(f  -  T0)  +  <p(t  -  T0)  +  tf0]. 
Using  the  trigonometric  identity 

cos  A  cos  B  =  -  [cos(A  —  B)  +  cos(A  +  B)\ , 
z 

together  with  the  assumption  that  p>(t)  is  a  low-pass  process,  the  filtered  mixer 
output  is  given  by 


x(t)  =  |  cos  [a;0T0  +  <p(t)  -  <p(t  -  r0)] ,  (54) 

where  t  >  tq.  Since  (pit)  is  not  strictly  a  low-pass  process  this  equality  is  not 
exact.  However,  for  an  oscillator  with  good  spectral  purity  <p(t)  is  very  nearly 
a  low-pass  process  so  that  equation  (54)  is  a  good  approximation.  Proceeding 
as  before,  define  the  complex  signal 


m 


__  p2[woro+^(*)-^(t-T0)] 


(55) 


where  t  >  r0.  Using  definition  (22)  together  with  the  fact  that  J?$(r)  is  real,  it 
is  clear  that 

f  E[X(t)]=ei“°rR*(T) 

[£[**(*)]  =  e-^orJ?$(r). 


Therefore,  the  mean  of  the  process  x(t)  is  given  by 


E[x(t)]=-E 


X(i)+X*(t) 


=  -  COS  (w0  7*0  )  i?4>  (  To  ) . 


(56) 


Substituting  the  result  (35)  for  the  Wiener  process,  one  obtains 

E[x(t)]  =  ^cos(uoT0)e~aT0^2.  (57) 

Due  to  the  phase  noise  of  the  source  oscillator,  the  expected  D.C.  output  level 
(s/2)  cos(a;oTo)  is  reduced  in  magnitude  by  the  factor  e_"r°/2.  However,  for 
the  small  delays  usually  encountered  in  practice,  the  inequality  a  To  1  is 

well  satisfyied  so  that  e~ar° /2  ~  1.  In  this  case,  this  factor  is  neglidgible  and 
equation  (57)  reduces  to  the  usual  result  obtained  in  the  absence  of  any  phase 
noise.  For  large  time  delays  such  that  ar0  1,  the  two  mixer  inputs  become 
uncorrelated  and  E[x(t)]  — »  0. 
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The  rest  of  this  section  will  be  devoted  to  a  calculation  of  the  autocorrelation 
function  of  x(t).  Using  the  complex  signal  X(t),  this  autocorrelation  function 
may  be  written  in  the  form 


X(t2)+X*(t2) 


Rx(t2,h)  =  E[a:(t2)a:(ti)]  =  ~^E 

=  ^{£[Jf(f2)Jt((1)]  +B[X(i2)X*((1)] 

+  E[X‘{t2)X{h)\  +  £'[X-!(2)A'"lf,i]} 

It  is  convenient  to  define  the  shorthand  notation 


X(U)  +  V*(U) 


} 


(58) 


f  A <pi  =  (fiti)  -  (f(t !  -  To) 

\  X<P2  =  <P{t2)  ~  V{t2  ~  To). 

Thus,  the  first  and  second  terms  in  equation  (58)  may  be  written  in  the  form 

£[X(<2  )*(*!)]  =  ei2woT°E[ei^+A^],  (59) 

E[V(t2)X*(fi)]  =  E[e*  AV2-A<Pi)j_  (60) 

These  expectations  may  be  evaluated  as  follows:  assume  that  1 2  >  t\  >  tq. 
There  are  two  cases  to  consider  (figure  4):  case  1  where  ti  <  t2  —  To,  and 
case  2  where  t2  —  t0  <  ti.  In  case  1,  and  A^2  are  independent  and 

identically  distributed  Gaussian  random  variables  with  mean  zero  and  variance 
oltq .  Therefore,  using  the  known  characteristic  function  for  a  single  Gaussian 
random  variable,  one  obtains 

£[ei(AW+AW)]  _  E[eiA^]E[eiA^)  =  e~ar° ,  1 12  -h\>  r0.  (61) 

Similarly,  the  random  variable  (— A(pi)  is  also  Gaussian  with  mean  zero  and 
variance  a  7o  .  which  implies 

E[c*'( Av2-a^)j  _  E[eiAv*]E[e~iAv']  =  e~aT° ,  |t2  -  tx\  >  r0.  (62) 

In  case  2,  Ay>i  and  A^>2  are  no  longer  independent.  In  this  case,  let 


{*1  =  <p(t 2  ~  Tq)  -  -  To), 

X2  =  -  <p(t2  -  To), 

Z3  =  ^(h)  - 
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Figure  4.  Illustration  of  overlap  between  the  phase 
increments  at  times  1  and  2 . 
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Because  ip(t)  is  a  Wiener  process,  the  quantities  aq,  X2,  and  £3  are  mutually 
independent  Gaussian  random  variables  with  mean  values  E[x\ ]  =  E[x2]  = 
E[xz\  —  0,  and  variances 


where 


=  a(r0  |r|)  •£  |r|  <  rQ? 

L  if  Irl  >  r0; 


sw  =  {  ^-w,  „ 

11  ii  n  >  to- 

Note  that  /(— r)  =  /(r)  and  g(r)  —  l//(r).  Substituting  equations  (67)  and 
(68)  into  equation  (58),  the  autocorrelation  function  of  the  filtered  mixer  output 
x(t)  is  found  to  be 

£2 

Rx{h,ti)  =  g-  e  ar°  [cos(2u;oTo)/(r)  +  £f(r)] .  (71) 

Since  this  function  only  depends  on  the  time  difference  r  =  f 2  — 1\ ,  the  process 
x(t)  is  wide-sense  stationary  as  required. 

Power  Spectrum  of  the  Mixer  Output 

Taking  the  Fourier  transform  of  the  autocorrelation  function  (71),  the  power 
spectrum  of  x(t)  is  given  by 


=  j  e~aT°  [cos(2lj0t0)F(u)  +  G(w)] ,  (72) 

where  F(uj)  and  G{u>)  are  the  Fourier  transforms  of  /(r)  and  g(r),  respectively. 
In  order  to  compute  these  transforms,  let 


/(r)  =  /(r)  -  1 

9{t)  =  9{t)  ~  1, 


so  that 


F(u)  =  F{u)  -  6{lo) 
G(u)  =  G(u>)  —  5(w). 


Using  equations  (69)  and  (70),  a  straightforward  evaluation  of  these  transforms 
yields 

p(u)  =  [cos(“To) “ iaTo>SJ^f- ~ e~"j  ■  <75> 


G(w) 


fcos(wro)  +  e+QT° 

a2+U2  [  (u;ro) 
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Note  that  G( ui)  is  obtained  from  F{uj)  by  simply  reversing  the  sign  of  a.  Sub¬ 
stituting  equations  (74),  (75),  and  (76)  into  equation  (72)  one  obtains 


S*(u)  =  j  e_ar° 


[cos(2a;oro)  +  l]<J(w)  + 


a2  +  u;2 


where 


A(u>)  =  cos(a;To)  [cos(2woro)  —  l] 

-  [cos(2w0ro)  + 1] 

+  e+ar°  —  e~ar°  cos(2woro). 


Using  the  identities 


ex  =  cosh(ar)  +  sinh(cr) 
i~x  =  cosh(:r)  —  sinh(:r), 


{1  —  cos(2x)  =  2  sin2  x 
1  +  cos(2x)  =  2  cos2  x, 

the  function  A(u)  can  be  written  in  the  more  compact  form 

A(u)  —  [cosh(oT0)  —  cos(iUT0 )]  2  sin2  [ujqTq  ) 

,  .  sin^ro)! 


+  sinh(o!ro)  —  (ar0)- 


(uT0) 


2  cos2(u?oTo). 


Discussion:  The  power  spectrum  (77)  is  the  sum  of  two  terms  where  the  first 
term  is  a  delta  function  and  the  second  term  is  a  continuous  function.  The  delta 
function  in  equation  (77)  cooresponds  to  the  power  contained  in  the  average 
D.C.  signal 

E[x(t)\  —  n  =  t  e~ar° /2  cos(w0t0).  (80) 


e  ar°/2  cos(uto'ro). 


For  the  purpose  of  studying  the  phase  noise,  the  delta  function  part  of  the 
power  spectrum  will  be  separated  from  the  continuous  part.  This  is  achieved 
by  separately  considering  the  power  spectra  of  the  two  signals  x(t )  —  jj,  and  \x  = 
const.  Mathematically,  this  amounts  to  using  the  autocovariance  of  the  process 
x(t)  instead  of  the  autocorrelation.  The  autocovariance  of  the  process  x(t)  is 
defined  by 

C(t2,ti)  =  [x(t2)  -  n(t2)]  [x(fi)  -  M(U)]  } 

=  R{h,h)  -  (81) 
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where  /j,(t )  =  E[x(t)\.  For  the  wide-sense  stationary  process  x(t),  the  mean  is 
constant  and  Cx(r)  =  Rx(t )  —  /i2  so  that  by  taking  a  Fourier  transform, 

Cx(u>)  =  Sx(oj )  -  n25(u).  (82) 

Hence,  the  delta  function  part  of  the  power  spectrum  has  been  subtracted  out 
and  Cx(ui)  is  simply  the  continuous  part  of  the  power  spectrum.  Obviously,  if 
the  process  mean  is  zero,  then  C(u>)  =  S(u).  For  the  mixer  output  x(t),  the 
continuous  part  of  the  power  spectrum  is  given,  using  equation  (77),  by 

(^)a(W).  (83) 

This  equation  is  one  of  the  main  results  of  this  report.  For  a  fixed  time  delay 
To,  the  first  two  factors  are  positive  constants.  The  frequency  dependence  of 
the  power  spectrum  is  contained  in  the  last  two  factors.  The  third  factor  is  the 
well  known  Lorentzian  spectrum  with  a  HWHM  given  by 

uub  =  <*  —  —,  (84) 

TC 

where  rc  is  the  correlation  time  of  the  oscillator  signal  u(t).  Comparing  this 
with  equation  (37),  the  HWHM  of  the  mixer  output  x(t)  is  twice  the  HWHM 
of  the  source  oscillator  v(t).  Therefore,  one  effect  of  the  mixing  process  is  a 
broadening  of  the  Lorentzian  part  of  the  power  spectrum  by  a  factor  of  2.  In 
addition,  and  more  importantly,  the  amplitude  of  the  spectrum  is  modulated 
by  the  factor  A{u>).  It  is  this  modulation  factor,  A(w),  which  is  responsible 
for  the  correlation  effect,  i.e.,  the  reduction  in  phase-noise  of  the  mixer  output. 
The  behavior  of  the  function  A(u)  is  investigated  in  appendix  B. 


Phase  Noise  of  the  Mixer  Output 

The  normalized  power  spectrum  of  the  mixer  output  x(t)  can  be  calculated 
from  the  PSD  Sx(oj).  Using  equation  (71),  the  total  average  power  in  the  signal 
x(t)  is  given  by 


Po  =  E[x2(t )]  =  Rx{ 0)  =  j  e~ar°  [e~aT°  cos(2 u>0t0)  +  e+ar°] .  (85) 

Applying  various  trigonometric  and  exponential  identities,  this  may  be  put  into 
the  form 


[cosh(aro)  cos2(u>ot0)  +  sinh(or0)  sin2(a>oTo)] . 


(86) 
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For  the  small  time  delays  met  in  practice,  it  is  almost  always  true  that  aro  <C  1. 
Therefore,  making  the  approximation  e+ar°  ~  e~aT°  ~  1,  equation  (85)  takes 
the  simplified  form 

e2  e2 

Po  ^  —  [cos(2cj0r0)  +  l]  =  —  cos2(a;oro).  (87) 

This  is  a  well  known  result  which  says  that  the  output  power  is  proportional 
to  cos2(u;o'ro),  essentially  the  cosine  squared  of  the  delay  (or  range).  Of  course, 
this  result  is  exact  when  the  phase  noise  of  the  source  is  zero.  But,  even  in 
the  presence  of  phase  noise,  this  result  is  a  very  good  approximation  as  long  as 
the  phase  noise  of  the  source  oscillator  is  low  and  the  time  delay  is  small,  i.e., 
ar0  <  1. 

In  order  to  calculate  the  phase-noise  spectrum,  the  delta  function  component 
of  the  PSD  (77)  must  be  omitted.  This  is  because  the  average  D.C.  component 
of  the  signal  x(t )  does  not  contribute  to  the  SSB  phase  noise  £(/).  Using  the 
continuous  part  of  the  PSD,  equation  (83),  together  with  equation  (12),  the 
normalized  phase-noise  spectrum  is  given  by 


2  Cx(u>) 
Po 


lu;=27r  / 


O?  +  U)‘ 


CD=27r/ 


where 


B  =  cosh(aro)  cos2(a>o7o)  +  sinh(aro)  sin2(u;oro).  (89) 

The  factor  of  two  on  the  right  side  of  equation  (88)  is  due  to  the  fact  that  both 
positive  and  negative  frequency  components  contribute  equally  to  the  measured 
noise  power.  The  function  A(u>)  may  be  written  in  the  compact  form 


AH  =  { [1  —  cos(wtq)]  +  a |  2  sin2 (uqTq) 


+  i  (aro)  1 


sin(wro) 

(WT0) 


+  b  >  2cos2(w0Tb), 


where 


{a  =  cosh(oro)  —  1 

b  =  sinh(aro)  —  (aro).  ^  ^ 

The  same  comments  made  for  equation  (83)  are  applicable  to  equation  (88). 
To  reiterate,  equation  (88)  shows  that  the  phase-noise  spectrum  of  the  mixer 
output  is  the  product  of  two  terms.  The  first  term  is  just  the  Lorentzian  spectral 
function.  If  equation  (88)  is  rewritten  in  the  form 


£(/)  =  ( 


a2  +  cj2 


‘2  A{u>) 

B 


w=2nf 
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then  by  comparison  with  equation  (41)  it  is  clear  that  the  linewidth  of  the 
Lorentzian  term  is  two  times  the  linewidth  of  the  source  oscillator.  This  is  one 
effect  of  the  mixing  process.  A  more  important  effect  arises  from  the  second 
term  2 A(u)/B  which  imposes  a  modulation  onto  the  Lorenzian.  It  is  this 
modulation  factor  which  is  responsible  for  the  coherence  effect,  i.e.,  phase- noise 
reduction. 


Phase  Noise  Reduction  Factor 

As  stated  in  the  introduction,  the  main  objective  is  to  compare  the  phase 
noise  of  the  mixer  output  in  figure  la  to  the  phase  noise  of  the  source.  The 
ratio  of  the  power  spectral  density  of  the  filtered  mixer  output  x(t)  to  that  of 
the  source  oscillator  y(t )  is  given  by 


K{u) 


sx(v) 

Sy  M' 


This  will  be  called  the  reduction  factor.  Substituting  the  power  spectra  (49) 
and  (83)  into  (93),  and  assuming  there  is  no  attenuation  of  the  delayed  signal 
(a  =  1),  the  reduction  factor  is  found  to  be 


K(u)  =  e- 


where  A(lo')  is  given  by  equation  (7.6).  This  equation  is  the  main  result  of  this 
report.  With  minor  modifications,  equation  (8.2)  may  be  written  in  the  form 


(aro/2)2  +  (urof 

(ar0)2  +  (coto)2 


2A(uj). 


For  the  time  delays  encountered  in  practice,  it  is  generally  true  that  ctTo  1. 
Consequently,  the  first  and  second  factors  in  equation  (95)  are  approximately 
unity,  and  it  is  a  good  approximation  to  write 

K(u)  ~  2 A{lo).  (96) 


Therefore,  the  frequency  response  of  K(u>)  is  contained  in  the  function  A(co) 
which  again  is  given  by: 

A(lo)  =  |  [l  —  cos(wto)]  +  aj  2  sin2 (wo to) 

+  |(ar0)  1 - (^V'f)]  +  6 1  2cos2(a;oro).  (97) 
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The  general  behavior  of  this  function  is  studied  in  appendix  B.  The  first  thing 
to  note  is  that  A(u>)  is  a  rapidly  varying  function  of  the  delay  r0.  Due  to  the 
factors  sin2(a»o7'o)  and  cos2(u;oro),  the  first  and  second  terms  in  equation  (97) 
oscillate  in  and  out  of  phase  as  To  cycles  through  a  period  of  l//o-  For  example, 
if  /o  =  33  GHz,  then  l//o  =  30  ps  (picoseconds)  which  shows  that  these  terms 
axe  oscillating  very  rapidly  indeed.  Therefore,  it  is  necessary  to  consider  three 
separate  cases: 

Case  (1):  u>oTo  =  mr,  n  =  1,2,3,...; 

Case  (2):  lo0t0  —  (2n  +  1)(7t/2),  n  —  0,1,2,...; 

Case  (3):  ujoTo  —  (2n  +  l)(7r/4),  n  =  0,1,2, _ 

Equivalently,  these  cases  are  characterized  by 

Case  (1):  cos2(woro)  =  1, 

Case  (2):  sin2(<x?oTo)  =  1, 

Case  (3):  sin2(o>oTo)  =  cos2(u;oTo)  =  1/2. 

The  behavior  of  A{u)  in  each  of  these  special  cases  is  discussed  in  appendix  B. 

To  illustrate  the  above  results,  the  reduction  factor  was  calculated  for  a  36 
GHz  oscillator  with  a  linewidth  parameter  r]  =  10-7.  For  the  purpose  of  making 
calculations,  it  is  convenient  to  introduce  the  variables 


X  =  (XTq 

y  =  wr0. 


Making  these  substitutions,  the  reduction  factor  (95)  takes  the  form 


2^). 

xz  +  yl 


where 


A(u)  =  |  [l  -  cos(t/)]  +  a |  2sin2(w0To) 


+  U*) 


1  _  sin (y) 


+  b  >  2cos  (wqto), 


(100) 


a  =  cosh(x)  —  1 
b  =  sinh(x)  —  x. 


(101) 


24 


In  the  low  frequency  limit,  as  w  -4  0,  equation  (107)  shows  that  K(u> )  -4 
(aro)2/2  where  the  additional  factor  of  1/4  comes  from  the  second  term  in 
equation  (94).  To  locate  the  frequency  of  transition  from  region  1  to  region  2, 
expand  the  cosine  in  equation  (107)  to  obtain 

K(lo)  «  2[(cjto)2  +  (aro)2],  (108) 

which  is  valid  for  lotq  small.  The  transition  occurs  when  utq  ~  aro  or  u  —  a. 
To  obtain  the  frequency  of  transition  from  region  2  to  region  3,  note  that  A(u>) 
reaches  its  maximum  at  loto  —  7 r,  i.e.,  u>  =  ~ / r0 .  where  the  maximum  value  is 
A(ujto  —  7r)  =  8  +  2(aro)2.  These  results  may  be  summarized  as  follows. 

Region  1:  0  <  u>  <  a,  K(u>)  ~  (aro)2/2  =  Kq ; 

Region  2:  a  <  u  <  7r/r0,  K{u)  ~  if0  +  2(u;r0)2; 

Region  3:  u>  >  7t/to,  RT(ci;)  ~  4. 

The  value  of  4  in  region  3  is  one-half  the  maximum  value  of  8. 

Also  plotted  in  figure  5  is  the  approximate  formula 

K(lo)  ~  2 [l  —  cos(a;ro)]  =  4sin2(u;ro/2)  (109) 

derived  from  FM  modulation  theory  (refs  1  through  3).  A  different  and  much 
simpler  derivation  of  equation  (109)  is  given  in  appendix  C.  It  is  clear  from 
figure  5  that  equation  (109)  is  a  good  approximation  in  regions  2  and  3  where 
it  approaches  the  exact  solution  for  case  3.  However,  (109)  fails  in  the  low 
frequency  limit,  that  is,  in  region  1.  The  breakdown  of  this  formula  can  be 
remedied  by  simply  adding  a  constant.  Thus,  based  on  the  previous  analysis,  a 
more  accurate  formula  is 


K(w)  £3  4[l  -  cos(u;ro)]  +  (ar0)2/ 2,  (110) 

where  a  =  dB  —  2r/ujQ  and  io^dB  is  the  HWHM  of  the  source  oscillator.  This 
formula  represents  a  slight  improvement  over  equation  (109). 

So  far,  the  discussion  of  the  solution  has  focused  on  the  three  special  cases: 
cases  1,  2,  and  3.  With  this  done,  it  is  necessary  to  discuss  what  happens  when 
the  delay  tq  lies  somewhere  in  the  range  between  these  cases.  It  is  possible  to 
write  equation  (97)  in  the  form 


A(lo)  =  A\  cos2(u>0t0)  +  A2  sin2 (ojqto), 


where 


A\ 


Aa(w)  =  2 


(ar0) 


sin(u;To) 

(wr0) 


(111) 

(112) 
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and 


A2  =  A2(u>)  =  2  j  [l  —  cos(o;to)]  +  a  j.  (113) 

The  function  A(ua)  reduces  to  A\(u)  in  case  1  and  A2{u)  in  case  2.  By  virtue 
of  the  identity 

sin2(0)  +  cos2(0)  =  1,  (114) 

equation  (111)  shows  that  the  response  A(u)  lies  between  A-y{ui)  and  A2 (w)  for 
all  values  of  tq  between  case  1  and  2. 


There  is  another  prominent  feature  of  figure  5,  namely,  that  the  separation 
between  the  response  curves  for  cases  1  and  2  is  independent  of  frequency.  This 
is,  in  fact,  a  general  feature  of  the  solution.  Ignoring  the  constants  a  and  b  in 
equations  (112)  and  (113),  the  ratio  of  the  response  in  case  1  to  that  in  case  2 
is  given  by 


M(u) 
A2  (w) 


~  ( ar0 ) 


1  -  sin(y)/y 
l-cos(y) 


(115) 


It  is  interesting  that  the  quotient  appearing  in  equation  (115)  is  very  nearly 
constant  on  the  interval  from  y  =  0  to  y  =  7r.  While  this  fact  is  difficult  to 
demonstrate  analytically,  it  is  easy  to  show  graphically  in  figure  6  where  the 
value  of  this  quotient  is  seen  to  be  approximately  1/3.  Therefore,  the  ratio  of 
Ai(u>)  to  A2(w)  is  approximately  constant  and 


Ai(w )  ^  cxtq 

a2(u)  ~  IT’ 


(116) 


where  ar0  =  2r)u>oTQ.  Thus,  for  the  delay  t0  =  10  9  seconds  in  figure  5,  the 
separation  between  the  response  curves  for  cases  1  and  2  is  aro/3  =  —48.2  dB. 

As  discussed  in  connection  with  equation  (97),  the  reduction  factor  is  a  very 
sensitive  function  of  the  delay  r0.  More  specifically,  as  To  cycles  through  a 
period  of  l//o  =  28  picoseconds  the  reduction  factor  goes  from  case  1  to  case  2 
and  then  back  to  1  again.  Of  course,  in  the  middle  is  case  3.  By  inspection  of 
figure  5  the  difference  between  case  2  and  case  3  is  3  dB.  Analytically,  in  case 
3  one  has,  using  equations  (111)  and  (116), 

Mu)  =  g  [M(v)  +  A2{u>)]  ~  -A2(u>),  (117) 

since  aro  -C  1.  This  accounts  for  the  3  dB  difference.  Note  that  the  response 
lies  between  cases  2  and  3  for  one-half  of  the  period  l//o  due  to  equation  (111), 
i.e.,  case  3  occurs  when  0  —  tt/2,  case  2  occurs  when  6  =  tt,  and  case  3  occurs 
again  when  9  =  3 tt/2.  The  question  arises:  for  what  percent  of  the  period  l//0 
is  the  response  within  3  dB  of  case  3?  The  answer  is  67  %  of  the  time.  To  see 
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Figure  6.  The  quotient  ( numer/denom)  in  equation 
(115)  is  approximately  constant  from  zero  to  pi. 


this,  note  that  if  6  —  0  corresponds  to  case  1  and  9  =  tt/2  corresponds  to  case 
3,  then  when  6  —  7r/3, 


A(u)  =  ^Ai(w)  +  ^A2(u)  ~  ^A2( u),  (118) 

which  corresponds  to  case  3  minus  3  dB.  Hence,  the  reduction  factor  is  within 
3  dB  of  the  solution  in  case  3  throughout  67  %  of  the  period  l//o- 

In  practice,  these  rapid  oscillations  of  the  reduction  factor  can  be  ignored 
since  in  a  radar  system,  for  example,  the  radar  return  is  “smeared  out”  in  time. 
That  is,  the  return  signal  is  the  sum  of  many  components  with  different  delays. 
This  results  in  an  averaging  effect  whereby  the  reduction  factor  is  averaged  over 
one  period  of  the  variable  0  =  uoTq.  Since  the  average  of  sin2(0)  and  cos2  (6) 
are  both  equal  to  1  /2,  the  averaging  effect  yields 

Mu)  =  \  {AAu)  +A2(uj)]  ~  ij42(w),  (119) 

Therefore,  using  equation  (96),  the  reduction  factor  in  the  presence  of  the  av¬ 
eraging  effect  is  given  by 

K(lo)  ~  A2(oj)  =  2[l  —  cos(u>r0)]  +  (cn-0)2,  (120) 

or,  equivalently, 

K(u) )  ~  4sin2(wro/2)  +  (aro)2.  (121) 


Conclusions 

The  phase  noise  of  an  amplitude  stabilized  oscillator  was  modeled  as  a 
sinewave  with  a  random- walk  phase  function.  The  propagation  of  phase  noise 
through  a  mixing  circuit  was  studied  to  determine  the  reduction  of  phase  noise 
due  to  mixing  the  time  delayed  signal  with  itself.  The  normalized  (baseband) 
power  spectrum  of  the  source  oscillator  is  given  in  equation  (41)  and  the  spec¬ 
trum  of  the  filtered  mixer  output  is  given  in  equation  (83)  or  (88).  The  main 
result  of  the  report  is  the  reduction  factor  (correlation  factor)  which  is  given  in 
equation  (95).  The  results  show  that  the  well  known  formula  (1)  for  the  reduc¬ 
tion  factor  is  a  good  approximation  except  in  the  low  frequency  limit  (o>  — >  0). 
An  improved  formula  was  given,  either  equation  (110)  or  (121),  which  is  valid 
for  any  frequency  offset. 

It  is  important  to  realize  that  the  theoretical  results  obtained  in  this  study 
are  only  applicable  to  the  extent  that  the  given  mathematical  model  is  valid.  In 
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general,  an  oscillator  will  exhibit  both  AM  and  FM  noise  components;  whereas, 
in  this  study  it  was  assumed  that  the  amplitude  noise  is  zero.  Nevertheless, 
for  an  amplitude  stabilized  oscillator  employing  current  state-of-the-art  tech¬ 
nology,  the  AM  noise  is  usually  negligible  compared  to  the  FM  noise.  Another 
limitation  of  the  model  employed  here  is  its  restriction  to  a  random  walk  phase 
function.  This  is  not  always  the  case  for  electrical  oscillators.  A  microwave  os¬ 
cillator,  for  example,  does  not  exhibit  a  perfect  random- walk  ( /“ 2)  phase  noise 
spectrum.  In  practice,  real  oscillators  have  phase  noise  spectra  which  exhibit 
combinations  of  different  slopes,  i.e.,  /~4,  /~3,  f~2,  /_1,  /°,  and,  more  gener¬ 
ally,  noninteger  slopes  f~a  where  a  >  0  is  a  real  parameter.  Nevertheless,  the 
random  walk  /~2  case  is  one  of  the  most  important  since  it  is  mathematically 
tractable  and  since  it  is  characteristic  of  a  resonantly  tuned  oscillator  operating 
near  its  thermal  performance  limit  (ref  12). 
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Appendix  A 

Necessary  and  Sufficient  Conditions 
for  Stationarity 
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For  a  zero  mean  Gaussian  phase  modulation  <p(t)  (not  necessarily  stationary) 
it  is  easy  to  obtain  a  necessary  and  sufficient  condition  for  the  stationarity  of 
the  process  v(t).  Since  the  sum  and  difference  of  Gaussian  random  variables  is 
Gaussian  the  increment  <p(t2)  —  p{t\)  is  a  Gaussian  random  variable  with  mean 
zero  and  variance 

v2(h,ti)  =  Var  | V(t2)  -  ^(*1)] 

=  E([p(t2)  -  ¥?(fi)]2)  =  cr%{t2)  -  2Rip(t2,t1)  +  cr^(ii). 

(A-l) 


By  using  the  characteristic  function  for  a  Gaussian  random  variable,  it  follows 
that 

R*{t  2,h)  =E(ei^~^t^  =  exp  [— cr2  (t2 ,  )/2] .  (A-2) 

Note  that  this  is  a  real  function.  Consequently,  using  equation  (24)  the  auto¬ 
correlation  function  of  v(t)  takes  the  form 


Rv(t  2,ii)  =  -jR^(t2,ti)cos(u;o^)- 


(A-3) 


Since  £J[u(f)]  =  0,  this  equation  implies  that  v(t)  is  wide-sense  stationary  if  and 
only  if  the  function  i?$(i2,ti)  depends  only  on  the  time  difference  r  =  <2  —  <1. 
From  the  previous  considerations,  it  follows  that  if  a2(t2,  h)  depends  only  on  r, 
then  i?$(t2,fi)  =  i2$(r)  is  strictly  a  function  of  r.  Conversely,  if  R$(t2,ti)  = 
R<s>(t)  is  strictly  a  function  of  r,  then  the  variance  cr2(t2,t  1)  is  also  a  function 
of  t.  This  proves  the  following  theorem. 

Theorem.  If  ip(t)  is  a  zero  mean  Gaussian  process ,  then  the  stochastic  process 
v(t)  =  cosfwot  +  pit)  +  <p 0]  is  wide-sense  stationary  if  and  only  if  Var\p{t2)  — 
p(ti)]  is  a  function  of  r  only,  that  is: 


Var[p(t2)  -  <p(h)]  =<jl(t2)-2Rv(t2,t1)  +  al(ti)  =  a2(T).  (A-4) 

In  this  case  the  autocorrelation  function  of  v(t)  is  given  by 

Rv(t)  =  ^i?$(r)cos(w0r),  (A-5) 

where 

R<s,{t)  =  exp[— i<72(t)].  (A-6) 

In  the  special  case  when  <p(t)  is  a  stationary  process,  (A-l)  is  only  a  function 
of  r.  Therefore,  one  has  the  following  corollary. 

Corollary,  If  <p(t)  is  a  stationary  zero  mean  Gaussian  process,  then  the  sto¬ 
chastic  process  v(t)  is  wide-sense  stationary  with  autocorrelation  function  given 
by  (A-5)  and  (A-6). 
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Appendix  B 

Analysis  of  the  Amplitude  Factor 


39 


The  functional  dependence  of  the  reduction  factor  K(u>)  is  primarily  con¬ 
tained  in  the  function  A(oj).  Therefore,  a  study  of  A(ui)  will  shed  light  on  the 
general  behavior  of  K(u>).  Since  A(u)  also  depends  on  the  time  delay  t0  it  will 
sometimes  be  denoted  by  A(u.  To),  where 

A(u>,to)  =  [cosh(ar0)  —  cos(wro)]2sin2(a;oT'o) 

+  sin 

It  is  clear  that  A( w,  To)  >  0  for  all  tq  >  0  and  for  all  u  real.  Also,  A{— u)  = 
A( oj),  so  that  A{uj)  is  an  even  function  of  lo.  One  important  observation  should 
be  made  at  this  point:  Due  to  the  factors  sin2 (wo^o)  and  cos2(woTo)  the  relative 
contributions  of  the  first  and  second  terms  alternate  as  To  cycles  through  a 
period  of  length  To  =  l//o-  For  frequencies  in  the  microwave  range,  To  is  on 
the  order  of  100  picoseconds  and  therefore  A(u)  is  a  very  sensitive  function  of  To- 
That  is,  the  first  and  second  terms  are  very  rapidly  going  in  and  out  of  phase. 
On  the  other  hand,  in  a  radar  system,  for  example,  the  radar  return  is  “smeared 
out”  in  time  so  that  on  average  each  of  these  terms  will  contribute  equally  to 
the  mixer  output.  This  averaging  effect  is  discussed  further  in  section  8.  In 
equation  (B-l),  as  To  increases  from  0+  there  are  three  cases  to  consider: 

Case  (1):  u>oTo  =  mr,  n  =  1,2,3,...; 

Case  (2):  uqtq  =  (2n  -f  1)(7t/2),  n  =  0,1,2,...; 

Case  (3):  u;0To  =  (2n  +  l)(7r/4),  n  =  0,1,2, - 

Equivalently,  these  cases  are  characterized  by 

Case  (1):  cos2(cjoTo )  =  1, 

Case  (2):  sin2^^)  =  1, 

Case  (3):  sin2(u;oTo)  =  cos2(u>oTb)  =  1/2. 

Because  the  mixer  output  x(t)  is  proportional  to  cos2 (woTo),  case  (1)  might 
be  expected  to  have  the  dominant  effect  on  the  output  phase  noise.  However, 
it  will  be  shown  that  the  coefficient  of  the  sin2  term  in  (B-l)  is  on  the  order 
of  unity,  while  the  coefficient  of  the  cos2  term  is  on  the  order  of  ar0  where,  in 
general,  aro  <§C  1.  Therefore,  cases  (2)  and  (3)  are  also  important. 

CASE  1:  wo^o  =  nn. 

In  this  case,  A(lo)  is  given  by 

A(uj)  =  2  sinh(aTo)  -  (ar0) - 1  °)  .  (B-2) 

\UTo)  . 


h(ar0)  -  (ar0) 


sin(u;To) 

(t*>r0) 


2  cos2  (wqTo). 


(B-l) 
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For  the  delays  usually  encountered  in  practice,  it  is  generally  true  that  ar0  1. 
Making  the  approximation  sinh(;r)  ~  x  +  rc3 / 3!,  this  takes  the  form 


A(to)  ~  2(ar0) 

where  b  =  (aro)3/3!.  The  magnitude  of  this  term  is  determined  by  the  factor 
ar0.  This  function  is  sketched  for  the  case  n  =  1  in  Figure  B-l.  Note  the 
linear  scale.  It  is  apparent  that  the  function  A(to)  will  cause  a  reduction  or 
nulling  out  of  the  phase  noise  spectrum  for  frequencies  in  the  neighborhood  of 
zero.  For  frequencies  from  zero  to  tor0  =  i r,  A(oo)  increases  monotonically  from 
the  minimum  value  2b  to  the  value  2 faro  +  b).  For  frequencies  greater  than 
uito  =  7T,  A{to)  is  nearly  flat  as  it  approaches  the  asymptotic  value  2(aro  +  b). 
Consequently,  the  nulling  effect  is  greatest  from  to  =  0  to  cotq  —  tt.  To  gain  some 
feeling  for  the  overall  magnitude  of  the  nulling  effect  suppose  that  ar0  =  10-4. 
Then  2b  ~  10-16  and  2(ar0  +  b)  10-4.  The  nulling  effect  caused  by  A(to)  is 
quite  strong  in  case  1. 

Consider  now  the  behavior  of  A(uj,t0 )  as  a  function  of  r0.  For  a  fixed  time 
delay  tq  =  mr/too  (n  fixed),  the  zeros  of  sin(curo)  occur  at  toro  =  mix,  or 


sin(u;To) 
(wr0)  I 


+  26, 


(B-3) 


u>  m 

too  n 


(B-4) 


where  m  =  0, 1, 2, ... .  The  behavior  of  A(uj)  for  different  values  of  the  delay  n 
is  illustrated  in  Figure  B-2.  In  general,  A(uj)  increases  monotonically  from  the 
minimum  value  2b  at  u>  —  0  to  the  value  2(ar0  +  b)  at  uj  =  lo0  jn.  For  to  >  u>o/n, 
A(uo)  is  nearly  constant. 


CASE  2:  ujqTq  =  (2 n  +  l){n/2). 

In  this  case,  A(u>)  reduces  to  the  form 

A{to)  =  2[cosh(aro)  —  cos(wro)] .  (B-5) 

For  the  small  delays  encountered  in  applications  it  is  almost  always  true  that 
aro  1.  Therefore,  making  the  approximation  cosh(x)  ~  1  +  x2/2,  equation 
(B-5)  takes  the  form 


A(  u>)  =  2  [l  —  cos(a;ro)]  +  2a,  (B-6) 

where  a  —  (aro)2/2.  This  function  is  plotted  on  a  linear  scale  and  displayed 
in  Figure  B-3.  This  function  exibits  large  amplitude  swings  as  a  function  of  to. 
As  to  increases  from  zero  to  lot q  =  tt  the  function  A[lo)  increases  monotonically 
from  its  minimum  value  2a  to  its  maximum  value  2  +  2a.  For  example,  if 
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Figure  B-l.  Linear  plot  of  the  amplitude  response  versus 
frequency  in  case  1. 


Figure  B-2 .  Amplitude  versus  frequency  in  case 
for  different  values  of  the  delay  n. 
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cxtq  =  10-4,  then  2 a  =  10-8  and  2  +  2 a  «  2.0.  Therefore,  in  this  example, 
A(uj)  spans  eight  orders  of  magnitude  between  to  =  0  and  ojtq  =  n.  Since  the 
phase  noise  spectrum  is  located  in  the  neighborhood  of  u>  =  0,  i.e.,  u  <C  u?o, 
this  shows  that  A(oj)  causes  a  significant  reduction  or  nulling  out  of  the  phase 
noise.  Of  course,  since  the  mixer  output  goes  as  cos2(u>oTo),  the  mixer  output 
is  zero  in  case  (2)  and  the  nulling  effect  is  irrelevant.  Nevertheless,  the  above 
analysis  is  relevant  for  a  study  of  case  (3). 

Next,  consider  the  dependence  of  A(u,  To)  on  the  time  delay.  In  case  (2),  the 
time  delay  is  given  by  r0  =  (2n  +  1)7t/2u;o.  The  extrema  of  A(lo)  occur  when 
ojto  =  miT ,  m  =  1,2,....  For  a  fixed  value  of  to  (fixed  n),  the  extrema  occur 
for  those  frequencies  which  satisfy 

oj  2m 
lo  o  2n  +  1  ’ 

where  m  =  1, 2, . . . ,  n  fixed.  The  behavior  of  A{ui)  for  different  values  of  n  are 
shown  in  Figure  B-4.  As  the  time  delay  n  increases,  the  bandwidth  from  ui  =  0 
to  the  first  maximum  u>/u> o  =  2/(2 n  +  1)  decreases.  Consequently,  the  band¬ 
width  of  the  nulling  effect  also  decreases  as  shown  in  Figure  B-4.  Nonetheless, 
the  nulling  effect  is  significant  except  for  very  large  values  of  n  (large  delays) 
since  the  phase  noise  is  concentrated  in  a  narrow  frequency  band  close  to  the 
carrier  frequency. 

CASE  3:  uqTo  =  (2n  +  l)7r/4. 

This  case  is  a  combination  of  cases  1  and  2.  Making  the  same  approximations 
for  sinh  and  cosh  as  in  equations  (B-3)  and  (B-6),  one  obtains 

A(u>)  —  (ar0)  1  -  +  [l  -  cos(wro)]  +  (a  +  b),  (B-8) 

where 

a  =  (ar0)2/  2,  6  =  (ar0)3/ 3!.  (B-9) 

Since  ar0  <C  1,  the  first  term  in  (B-8)  is  much  less  than  the  second  term  and 
the  overall  response  is  very  similar  to  that  of  case  2  shown  in  Figures  B-3  and 
B-4. 
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Appendix  C 

Derivation  of  the  Approximate  Formula 
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In  this  appendix  an  elementary  derivation  is  given  of  the  approximate  formula 
(1).  The  phase  of  the  output  signal  is  (f>{t)  —  <p(t )  —  <p(t  —  To),  where  p{t )  is 
the  phase  of  the  source.  Therefore,  the  autocorrelation  function  of  the  output 
phase  is  given  by 

S4>(r)  =  E[4>(t2  )4>(t1)] 

=  #{  M*2)  -  ^(*2  -  To)]  -  To)]  } 

=  2  Rv{t)  —  Rv{t  +  To)  —  Rv(t  —  To).  (C-l) 

Taking  the  Fourier  transform  and  using  the  time  shift  property,  this  yields  the 
result 

S^(w)  =  2  [l  -  cos(a;ro)]5v,(a;).  (C-2) 

Therefore,  the  ratio  of  the  power  spectrum  of  the  output  phase  6{t)  to  the 
power  spectrum  of  the  input  phase  p{t)  is  given  by 

7-^7— [  =  2[l  -  cos(a;ro)]  =  4sin2(o;ro/2).  (C-3) 

This  is  just  formula  (1).  The  problem  with  this  result  is  that  the  power  spectrum 
of  the  oscillator  signal  u(f)  is  not  equal  to  the  power  spectrum  of  the  phase 
function  pit),  and  the  PSD  of  x(t)  is  not  equal  to  the  PSD  of  <f>(t).  In  the  case 
of  a  random-walk  phase  function,  a  more  careful  analysis  leads  to  the  correct 
result  (94). 
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